home *** CD-ROM | disk | FTP | other *** search
/ NeXTSTEP 3.3 (Developer)…68k, x86, SPARC, PA-RISC] / NeXTSTEP 3.3 Dev Intel.iso / NextDeveloper / Headers / g++ / HypGeom.h < prev    next >
C/C++ Source or Header  |  1995-02-06  |  2KB  |  72 lines
  1. // This may look like C code, but it is really -*- C++ -*-
  2. /* 
  3. Copyright (C) 1988 Free Software Foundation
  4.     written by Dirk Grunwald (grunwald@cs.uiuc.edu)
  5.  
  6. This file is part of the GNU C++ Library.  This library is free
  7. software; you can redistribute it and/or modify it under the terms of
  8. the GNU Library General Public License as published by the Free
  9. Software Foundation; either version 2 of the License, or (at your
  10. option) any later version.  This library is distributed in the hope
  11. that it will be useful, but WITHOUT ANY WARRANTY; without even the
  12. implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
  13. PURPOSE.  See the GNU Library General Public License for more details.
  14. You should have received a copy of the GNU Library General Public
  15. License along with this library; if not, write to the Free Software
  16. Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.
  17. */
  18. #ifndef _HyperGeometric_h
  19. #ifdef __GNUG__
  20. #pragma interface
  21. #pragma cplusplus
  22. #endif
  23. #define _HyperGeometric_h 
  24.  
  25. #include <Random.h>
  26.  
  27. class HyperGeometric: public Random {
  28. protected:
  29.     double pMean;
  30.     double pVariance;
  31.     double pP;
  32.     void setState();
  33.  
  34. public:
  35.     HyperGeometric(double mean, double variance, RNG *gen);
  36.  
  37.     double mean();
  38.     double mean(double x);
  39.     double variance();
  40.     double variance(double x);
  41.  
  42.     virtual double operator()();
  43. };
  44.  
  45.  
  46. inline void HyperGeometric::setState() {
  47.   double z = pVariance / (pMean * pMean);
  48.   pP = 0.5 * (1.0 - sqrt((z - 1.0) / ( z + 1.0 )));
  49. }
  50.  
  51. inline HyperGeometric::HyperGeometric(double mean, double variance, RNG *gen)
  52. : Random(gen) {
  53.   pMean = mean; pVariance = variance;
  54.   setState();
  55. }
  56.  
  57. inline double HyperGeometric::mean() { return pMean; };
  58.  
  59. inline double HyperGeometric::mean(double x) {
  60.   double t = pMean; pMean = x;
  61.   setState(); return t;
  62. }
  63.  
  64. inline double HyperGeometric::variance() { return pVariance; }
  65.  
  66. inline double HyperGeometric::variance(double x) {
  67.   double t = pVariance; pVariance = x;
  68.   setState(); return t;
  69. }
  70.  
  71. #endif
  72.